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ANGULAR MOMENTUM DESATURATION FOR SKYLAB 
USING GRAVITY GRADIENT TORQUES 

SUMMARY 


The angular momentum desaturation method for the Skylajb is presented. This 
method utilizes the gravity gradient torques and therefore minimizes the necessity for mass 
expulsion by the thruster attitude control system (TACS). The desaturation method 
requires maneuvers about the two principal axes of large inertia. The percentage of the 
orbit used for desaturation is selectable. An arbitrary misalignment between the axes of 
principal moments of inertia and the geometric vehicle axes is permissible. An angle bias 
about the sun line minimizes the momentum accumulation in the orbital plane. 

This report is an extensive revision of Reference 1 . The desaturation scheme has 
been expanded to include arbitrary principal moment-of-inertia axes misalignment, adaptive 
maneuver limiting according to the available maneuver momentum, and reduction of the 
third-order sampled data system to a second order. 


INTRODUCTION 


The Skylab Apollo Telescope Mount (ATM) experiments require that the solar 
instruments remain inertially fixed (sun-oriented) during the day portion of the orbit. 

Gravity gradient, aerodynamic, venting, and other external torques acting on the vehicle 
during this time must be absorbed by an angular momentum storage device; in this instance, 
a system of three double-gimbaled control moment gyros (CMG’s) [ 2, 3, 4, 5] . Portions of the 
disturbance torques are noncyclic and tend to saturate the CMG system, which has a limited 
momentum storage capacity. A method for momentum desaturation that does not require 
mass expulsion is desired. The gravity gradient torques acting on Skylab are developed for 
small deviations from the sun-oriented reference coordinate system. These equations are 
used to show that maneuvers about the two axes of large moments of inertia are sufficient 
to desaturate the accumulated momenta about all axes. All attitude maneuvers for 
desaturation are made during the night portion of the orbit (unless an insufficient night 
portion is available, where part of the daylight portion is used), and the percentage of the 
orbit utilized for desaturation is selectable. Expressions for the desaturation angle 
commands are developed assuming infinite angular vehicle rates. The effect of finite rates 
then introduces an efficiency, which is calculated and compensated for by a change in the 
commanded maneuver. An arbitrary misalignment between the principal moment-of-inertia 



axes and the geometric axes of the vehicle is acceptable. An angle bias about the sun line 
allows the minimization of the momentum component along the projection of the Zp axis 

into the orbital plane (see Appendix A for coordinate system definitions, etc.). 

The desaturation angles have been assumed to be small enough to make small angle 
approximations valid which allows the principle of superposition and also allows the treat- 
ment of the angles as if they were vectors. 

To avoid endless repetition in the following discussion, “momentum” is used for 
“angular momentum,” “desaturation” is used for “angular momentum desaturation,” and 
“vehicle axes” means “geometric vehicle axes.” 


GRAVITY GRADIENT TORQUE 


The gravity gradient torque acting on the vehicle can be expressed as [6] 


l g = 3 a»Ir ] [I] r 


( 1 ) 


when the vehicle is in a circular orbit (which is assumed in the further development); [I] is 
the vehicle inertia matrix, f2 is the orbital rate, £ is a unit vector parallel to the radius 
vector from the earth center to the vehicle center of mass, and [r ] is defined as 



0 

‘ r z 

+r y 

[ 7 ] = 

+r z 

0 

- r x 


' r y 

+r x 
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When the torque is expressed in the principal reference system, PR (see Appendix 
A for definitions), when [Ip] is the principal moment-of-inertia matrix, and when [e] is 

the transformation from the PR system to the P system, the equation (1) becomes 


T g = 3fi 2 [rp r ][e] T [I p j [e] r pr 


( 2 ) 
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or 
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Aa = +c2r? (1 +c2r? t( j) 

A 12 = -si? s2rj td 
A i 3 +ct? s2r? td 

A 22 = -0.5 (1 + c2tj ) (1 + c2r? td ) + (1 - c2r? td ) 
A 23 = -0.5 s2r? (1 + c2r? td ) 

A 33 = +0.5 (1 - c2r?) (1 + c2i? td ) - (1 - c2r? td ) 


(3) 


(4) 


(5) 
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Appendix B gives a detailed derivation of equations (4) and (5) for the case where 
the Xp r axis is in the orbital plane, the Zp f axis has an elevation angle of i? from the 

orbital plane and is generally pointing toward the sun (the z y axis is pointing exactly 

toward the sun), 1 and the e angles are small. The gravity gradient torque has been split 
into a nominal part Tg n which is not a function of the e angles and a controllable part 

Tgd which will be used for the desaturation method. The nominal part shows that a bias 
momentum accumulates about the x pr axis (Fig. 1); the others have only cyclic terms 
(Fig. 2). Visualization of T gn as a function of time for an arbitrary moment-of-inertia 
distribution is described in Reference 7. 



MOMENTUM DESATURATION METHOD 


General 

The development of a desaturation method using the gravity gradient torques 
consists of maneuvering the vehicle through angles (e’s) in such a way that the momentum 
accumulation caused by the angles desaturates the stored momentum to keep the total 
momentum bounded and to avoid the need for desaturation by the thruster attitude con- 
trol system (TACS). The CMG attitude control system executes the e angle commands. 


1. r 7 t£ j is the orbital angle from desaturation midnight (to be explained later). 
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Figure 2. Typical cyclic angular momentum. 

This system is described in References 2, 3, 4, and 5. The desaturation maneuvers 
must be consistent with the mission constraints which for Skylab allow maneuvers for part 
or all of the night portion of the orbit only. During the daylight part of the orbit, the 
e angles must be zero except an angle about the vehicle z axis (sun line) which must be 
kept constant; therefore, the desaturation loop is only closed on a per orbit basis. The 
desaturation method can be separated into two parts. One part consists of the development 
of the desaturation momentum commands, using the total momentum profile. The other 
part consists of the generation of desaturation angles e, given the desaturation momentum 
commands. The latter affects the former and is therefore presented first. 


Desaturation Angle Generation 

The angular momentum desaturated by a given set of e angles is developed first. 
Since this momentum should be equal to the command, the equations are then inverted. 
Inspection of the Ay’s reveals that they are either even or odd functions of j?^. There- 
fore, assuming that we have one set of constant angles before midnight, e_ + Ae, another 
set after midnight, e_ - Ae, and the integration interval is from r? tc j = -p to r? tc j = +p, we 
get (Appendix C) 
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an = c2i?(2p+s2p) 

3l2 = ST? { 1 - c2p ) 
a i3 = -cr? (1 - c2p) 

a 2 2 = -0.5 (1 + c2t?) (2p + s2p) + (2 p - s2p ) 

a 23 — -0.5 s2r? (2p + s2p) 

a 33 = 0.5 (1 - c2t?) (2p + s2p) - (2p - s2p) 

Aj = 2/(3S2 AI X ) 

A 2 = 2/(3S2 Al y ) 

A 3 = 2/(3S2 AI z ) 

The moment-of-inertia distribution of the Skylab was configured so that Al x is small, 
resulting in a small bias momentum accumulation. But a small Al x makes the use of an 
e x or Ae x very ineffective for momentum desaturation, leading to the conclusion that no 
maneuvers about the Xp axis will be made, and we have as components of 
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H(j x — -ai 2 Aey/(-A2 ) + ai3 Ae z /A3 


(7) 



Equation (8) can be rewritten in the following form: 


( 8 ) 



Equation (9) shows that resolution into the orbital plane decouples the equation. 

It has been assumed so far that the e angles can be reached instantaneously. In the 
following, the angle profiles of Figure 3 are substituted. The difference between the ideal 
and the actual momentum desaturated is expressed by an efficiency, and we have (the 
components are in the orbital coordinate system CS Xp r ) 

(He, ■ E x (H dx ). d 

( Hdy )act = Ey ( Hdy )id 

(He ■ MH d 

where 


E x = (2sp - s2p)/[p(l - c2p)] 

E y = (cp - c2p)/(p s2p) 

E z = [3p/2-(cp -c2p)/p]/(2p -s2p) 

Appendix D gives a detailed development of the E’s. Equations (7) and (9) become 
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H dx = -I(2sp - s2p)/p] [stj Ae y /(-A 2 ) +ci? Ae z /A 3 ] 


( 10 ) 


H dy 


Cl? 

ST? 


a y 0 


C T? -ST? 


e y K- A a ) 

H d z 


-ST? 

CT? 


j 

C3 N 

O 

1 


-3 

O 

1 


e z^ 3 


with 


(ID 


a y = 2(cp - c2p)/p 
a z = -3p/2 + (cp - c2p)/p 


Before inverting equations ( 1 0) and (11) it is necessary to discuss the orbital desaturation 
parameters a y and a z . The behavior of ay for varying desaturation percentages (or 

varying p), normalized with respect to its maximum value, is shown as 


(cp - c2p)/p 
MAX [(cp -c2p)/p] 


in Figure 4. For a given orbital e y the maximum momentum is desaturated at a desatura- 
tion percentage of 32. 1 percent (p = 57.8°). Instead of assuming a given angle we can 
assume that a constant maneuver momentum is available and that the maneuver angle is 
proportional to the available desaturation time. This is shown as 


cp - c2p 

MAX (cp - c2p) 


The peak of the orbital y-momentum desaturation shifts to a desaturation percentage of 
42 percent (p = 75.5°), which is higher than the maximum nighttime available on Skylab 
(slightly less than 40 percent). 

A comparison of the orbital z-momentum parameter a z shows that even at the 

maximum desaturation percentage (about 40 percent) the same angle would only desaturate 
about one fourth in z as would in y. On the other hand, a constant angle about the 
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63 ° 



Figure 4. Relative desaturation capacity about orbital y axis. 


z axis (if held for one orbit) is more than three times as effective as the same angle in y. 
Since mission constraints allow a constant angle about the vehicle z axis (which desaturates 
only orbital z momentum), it will be used exclusively for orbital z desaturation (the 
appropriate relationships are developed later). 


Equations (10) and (11) show the actually desaturated momentum for a given set 
of desaturation angles and the maneuver profile of Figure 3. Assuming that the actual and 
the commanded momentum are equal, we can invert these equations. However, one more 
constraint must be added to equation (10). If we minimize [Aey/(-A 2 )] 2 + [Ae z /A 3 ] 2 
(minimum maneuver momentum for Iy = I z ), inversion yields 


Ae y /(-A a ) 

Ae z /A 3 


^dx 

2s p - s2 p 



( 12 ) 
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Inversion of equation (11) results in 


e y /(-A 2 ) 


C7J 

ST? 


l/*y 

1 

0 


CT? -ST? 


H dy 

€ Z /^ 3 


-ST? 

CT? 


0 

l K 


ST? CT? 


H dz 

- 


L J 



- 


_ - 




The assumption that all of the orbital z momentum, H dzo = ST?H dy + C7?H dz , will be 
desaturated by a rotation about the vehicle z axis allows us to set l/a z to zero; thus 


e y /(- A 2 ) 
e z /A 3 


t(ci?H dy - sr?H dz )/a y ] 


ci? 

-Sr? 


(14) 


The effectiveness (per orbit) of a rotation Av z about the vehicle z axis must be calculated 

to justify this assumption. Evaluation of equation (9) for p = n yields this effectiveness, 
where it is again assumed that Al x ' is negligibly small (no H dx ): 


‘ H dy 


CT? 

ST? 


'0 

0 ~ 


CT? 

-ST? 


’i/(-a 2 ) 0 ~ 


‘k 22 

K23" 


’ 0 " 

H dz 

Au z - 

-ST? 

CT? 


0 

-2-n 


ST? 

CT? 


0 i/a 3 


k 32 

K33 


Au z 



L J 


_ 

jm 


_ 



_ _ 



- 


. 


or 


H dzo = -27r [K 2 3 St?/(-A 2 ) + K 33 ct?/A 3 ] Av z (15) 

where 


[K] 


K„ K n K 13 

K21 K22 K23 

K31 K32 K33 


is the transformation matrix from vehicle axes to principal axes. Equation ( 15) shows that 
Au z about the vehicle z axis is still only effective for orbital z-momentum desaturation, 
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when considered on a per orbit basis, and does not couple into orbital y-momentum 
desaturation (which was to be shown). Some difficulties remain in the area of command 
generation which will be discussed later. Inversion of equation (15) results in 


Av z = -H dzo /|27r[K 2 3Srj/(-A 2 ) + K 33 CT 7 /A 2 ] 

= - (si?H dy + ChH dz )/ |27r [ K 23 stj/(-A 2 ) + 

The following approximation can be employed: 

cr? x /A 4 = K 23 sr?/(-A 2 ) + K 33 ct?/A 3 

which assumes that the moment-of-inertia misalignment is mostly about the x axis 
(Kji ~ 1 ) and that I y ~ I z such that an average can be used. 

A 4 = 0.5 [(-A 2 ) + A 3 ] (17) 

In summary, given three momentum components H dx , H dy , and H dz in principal 
coordinates, the desaturation angle commands in vehicle coordinates are 



Al) z - -A 4 (SJ?H dy + CTjH dz )/(27TC77 x ) 



Ku K 21 

k 12 k 22 

k 13 k 23 


k 31 " 


0 


k 32 


st?(-A 2 ) 

[-p/(2sp - s2 p)] H dx 

k 33 


cr?A 3 



e x 


‘ K„ K 21 K 31 ' 


0 

e y 

= 

K 12 k 22 k 32 


C7?(-A 2 ) 

e z 


K 13 K 23 K 33 


-stjA 3 

_ 


- 




p(cr)H dy - s??H dz ) 
2(c p - c2 p) 


(18) 


(39) 


( 20 ) 
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Equations (19) and (20) assume that the desaturation angle commands are applied accord- 
ing to the profiles in Figure 3 and equation (18) assumes that Au z is applied for a full 

orbit (where the change in u z is made right after the desaturation interval at *?td = » 

see Appendix E for nominal v z angle). 

The maneuver commands before midnight are then (these angles are reached at 

’ltd = ' p l 2) 


e xl 


‘•x 


Ae x 

e y 1 

= 

e y 

+ 

Ae y 

e zl 


e z 


Ae z 


The commands after midnight are (these angles are reached at 7}^ = p /2) 


( 21 ) 


e x2 


V 


Ae x 

e y2 

= 

e y 

- 

<1 

e z2 


e z 


Ae z 


( 22 ) 


The maneuver commands of equations (21) and (22) disregard the fact that the CMG’s 
may not have sufficient momentum reserve available to execute the maneuver which will 
require reduction of the maneuver angle commands. This is treated in Appendix F. 

Figure 5 shows a typical CMG momentum profile. 


Angular Momentum Sampling 

An adequate representation of the momentum accumulation of the Sky lab is 
H t = H s s2(r? t( j - 7r ) + H e c2(i7 t( j - it ) 

+ Hk(7?td- 77)/( 2 77) + (Ha + Hb) (23) 
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Figure 5. Typical CMG momentum profile. 
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where 7) t( j = J2t, H g and H c are the amplitudes of the cyclic components (see Appendix G), 
is the momentum per orbit caused by a constant torque or constant angle, and H a is 
the average momentum at noon in excess of the desired bias [equation (23) is set up 
with respect to noon and r? t( j is with respect to midnight] . 

The task of the desaturation method is to center the cyclic components about some 
given bias point (Hjj), which is normally zero. This necessitates a determination of H a 

and H^. To accomplish this, the total system momentum is sampled at four points in orbit 

as shown in Figure 6. The total system momentum is used rather than the CMG momentum 
only to avoid invalid readings in case the vehicle is maneuvering at the time the sample is 
taken. The sample points are at 2(7r/4), 3(n/4), 5 (tt/ 4), and 6 (tt/ 4) and accordingly 
carry the subscripts 2, 3, 5, and 6. 
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Evaluation of equation (23) at the sample points results in 


" — 12 " 


0 -1 -0.25 +l' 


'n s 

— 13 


- 1 0 -0.25 +1 


— c 

its 


+1 0 +0.25 +1 


& 

16 . 


0 -1 +0.25 +1 


- — a + — b. 


( 24 ) 


Inversion of equation (24) yields 


— s 


' +0.25 

-0.5 

+0.5 

- 0.25' 


12 " 

— c 


-0.5 

+0.5 

+0.5 

- 0.5 


— 13 

— k 

1 

-2 

0 

0 

+2 


— 15 

. Ha + H b - 


_ 0 

+0.5 

+0.5 

0 


H t6 - 


It should be recognized that as far as and (H a + H^) are concerned, only the spacing 
of the sample points is important; their relationship to r? tc j is immaterial. 


Desaturation Command Generation 

A desaturation command of the form 



(26) 


would result (in steady state) in an offset of H a = since must be desaturated each 
orbit. The following form would eliminate (for the ideal case) the offset: 


H d = - (H a + H k ) 


(27) 


The nonideal case (tolerances in the momentum measurements, etc.) still results in an offset. 
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A form that eliminates steady state offsets is 


Od = SHa-O k 


( 28 ) 


with 


2H a = 2H a(n ., ) -H a +0.75H a(n ., i 


(29) 


where (n-1) indicates the value for the past orbit. Any steady state signal needed will now 
be provided by 2H a , a computed quantity, rather than by a physical offset (H a ). The 

sampled data characteristics of equations (28) and (29) are treated in Appendix H. 

The fact that a u z -angle is utilized for desaturation of the orbital angular momentum 

component introduces a difficulty: A tamp in the orbital z momentum must be considered 
as having been introduced by the desaturation itself and its effect on H az must be com- 
pensated. A simple example will illustrate the point. Consider only an initial momentum 
about the orbital z axis. This momentum is measured as (H az ) j = H 0 during the first 

orbit. At the end of the first desaturation interval, u z is changed by Au z so that H 0 
is eliminated after exactly one orbit. But at noon of the second orbit, v z has not done its 
job yet, and (H az )2 measures the remaining momentum (about 0.55 H 0 for a 35 percent 
desaturation). This remaining momentum is 

<H rz > 2 = -0 55 (H kz >2 [since (H tz ) 2 = -(H az ),] 


and it must be subtracted from the orbital z momentum: 


(H az > 2 = y < H tz3 + H tz5>2 + 0 555 ( H kz>2 < 30 > 

For the example (H az )2 would be zero and (^2)2 would have the right value to eliminate 
exactly the previously introduced Av z . Equation (30) seemingly contradicts the comments 

made with respect to equation (26), but it should be recognized that in steady state there 
is no Hj coz ; i.e., whatever average torques are acting about orbital z are compensated for 

by the appropriate v z . 
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In a sense, the latter fact is a gravity gradient sensor holding the vehicle in a certain 
attitude about the sun line. It holds this attitude also if there is rate gyro drift, which is 
compensated for on a per orbit basis by the appropriate Av z . Since the external torques 

(besides gravity gradient) are unknown, it is impossible to state what attitude is maintained. 
For attitude control, a combination rate gyro/strapdown reference calculation is sufficient, 
but a star tracker or other reference is necessary for exact knowledge of the attitude about 
the sun line. 

The desaturation momentum commands used in equations (18), (19), and (20) are 
in principal coordinates, but the angular momentum samples and the desired momentum 
bias are in vehicle coordinates. Therefore, 


H k = 2[K] (H t6 - H t2 ) (31) 

H a = [K] [0.5 (H t3 + H t5 ) -H b ] + 0,555 [tj] [0,0,(sr?H ky + cr?H kz )] T (32) 
2H a = 2H a -H a + 0.75H a(n . 1) (33) 

H d =s Ha-H k (34) 


So far it has not been considered that the CMG’s can saturate, but if they do, a minimum 
impulse bit (MIB) must be fired in the appropriate direction by the TACS to desaturate 
the CMG’s immediately or fine pointing is lost. Therefore, the momentum desaturated by 
the MIB’s must be included in the samples, or H a and H k give erroneous results: 


H ti = H t -2H mib (i = 2, 3, 5, 6) 

where is the total system momentum at the sample time and 2H m j b is the accumu- 
lated MIB momentum fromthe time of the first sample (i = 2) to the time of the present 
sample. H k can then be calculated as usual [equation (25)] , but H a becomes 

H a = 0.5(H t 3+H t5 )-H b + SH mib (35) 

where £H m j b is the accumulated MIB momentum from the first sample until the time of 

execution of equation (35); thus whatever momentum the MIB’s have already desaturated 
does not need to be desaturated by maneuvers. 
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APPENDIX A 


DEFINITIONS OF SKYLAB COORDINATE SYSTEMS AND ANGLES 


Only the coordinate systems (CS’s) needed for the development of the desaturation 
method are defined. Some transformations are shown for completeness, but are not needed 
for the development, which is pointed out in the definitions (e.g., for X r : 77 y -angle is not 

needed). All CS’s are right handed and orthogonal. All angles are defined mathematically 
positive and those with a t subscript are orbital y angles. 


Symbols 

Transformation Matrix 

Definition 

x o 


Basic orbital CS. z Q toward ascending node; 

y 0 

N/A 

y 0 toward orbital north. (CS X Q is not 

z o 


needed explicitly for desaturation method.) 

x r 


Reference CS. z f toward sun; x r in the 

y r 

x r = [t? x ] [7? y ] x 0 

orbital plane; y r in northern orbital hemis- 

Z r 


phere. (17 y angle is not needed.) 

t-4 

O 

X 


Orbital reference CS. z Qr along projection 

y Q r 

X or " 1 * 1 x 1 T Xr 

of z r into orbital plane; y or toward 

z or 


orbital north. 

x d 


Disturbance CS. z^ toward center of earth; 

yd 

Xd — XqJ. 

y<j toward orbital north. (ij t = 0 indicates 

z d 


orbital midnight; is a y angle.) 

x vr 


Vehicle reference CS. For no attitude 

y,, r 

Xy r = I»zl X r 

deviation the vehicle geometric axes will be 

J Vl 

z vr 

aligned with this CS. 
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Symbol Transformation Matrix 


Definition 


Vehicle geometric CS, z y toward the ATM 
X v = [0 e ] X vr rack; x y toward CSM. [<p e ] is the attitude 

deviation (commanded or error.) 


P 

P 

P 


X p = [K] X y 


Principal axis CS. Axes along principal 
moment-of-inertia axes; labeling of the axes 
such that tr[K] is maximized. 




= [K] X w 


Principal axes reference CS. This is the CS 
for the principal axes in the absence of an 
attitude deviation. 



Orbital principal CS. z op along the pro- 
X op = [a zp ] T [r?] T X pr jection of z pr into the orbital plane; y op 

toward orbital north. (i> zp is not needed.) 


x pr 



Another definition of the CS X pr showing 

X pr = t^xp 1 f u zel l^tml X or the an § les ^tm and u ze which 310 Used in 

Appendix G. (p xp is not needed.) 
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APPENDIX B 


DERIVATION OF GRAVITY GRADIENT TORQUE 


The gravity gradient torque acting on the vehicle is [equation (2)] 

Ig = 3 ft 2 [r pr ] [e] T [I p ] [e]rp r (Bl) 

with 


0 

orbital rate for circular orbit 

unit vector parallel to radius vector from the center of the earth to the 
vehicle center of mass 

transformation matrix from the principal reference CS Xp r to the 
principal axes CS X p 

In the following development it is assumed that the e angles are small and that 
Xp r lies in the orbital plane and z pr points toward the center of the sun. Equation (2) 

can be written as 

Ig = 3ft 2 [r pr ] j [E] + [e] J [I p ] j [E] - [ 7 ]^ (B2) 


LipJ 


where 


ft 


r 

^P r 


[e] 


0 0 


I 


y 

o i 


o 


principal moment-of-inertia matrix 
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with 


m = 


0 

• e z 

+e 

y 

+e z 

0 

" e x 

' e y 

+e x 

0 


(B3) 


where e x , Cy, and e z are small angles about the corresponding principal axes and [E] 
is the identity matrix. 

The moment-of-inertia differences will appear so frequently in the following that 
new symbols are introduced: 


A *x V* *y 
A *y ~ ^x ' ^z 


Al z ly * I x 


Operating on equation (B2) leads to 


Tg = 3ft 2 [r pr ] [I p ] - [I p ] [e] + [T] [I p ] - [7] [I p ] f?] _r pr 


The term [e] [Ip] [7] contains only square terms in e and is neglected. 


T„ = 3ft 2 


0 -r„ 


+r„ 


+r„ 


+r. 


0 


+I x 

-Al z ey 

i 

t-H 

~^z e z 

+I y 

’ Al x e 

-AlyEy 

- Al x e x 

+I Z 
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where r x , r y , and r z are the components of rp r in the principal reference CS Xp r 


Al x r y r z 


" r y +r z ' r x r y +r x r z Al x e x 


T = 3J2 2 / AI y r z r x + +r y r x -r|+r x 2 -r y r z 


Aly e y 


A *z r x r y 


+r r 
Vy 


-r 2 -f-r 2 AT e 

1 X Tiy ^1 Z e z 


The torque has been split into two parts: one that is independent of the e angles and can 
be considered the nominal part, and one that is dependent on the e angles and can be 
controlled. 

The components of rp r are developed as follows: 


Ipr = t 7 ?] ^tm 1 l^t^id = [7?1 lr? tdl T Id 


[r? td ] = 


r x 1 0 0 C7 ltd 0 S7 ?td 0 


r y = 0 ct? st? 


0 1 0 


0 -sr? ct? _ST ?td ^ CT ? td '1 


S7 ?td 

ST? CT? td 
CT? CT? t( j 
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Consequently, 


T = 
-g 


Xgn + Igd 


Ign = 30* 


AI x st? ctj c 2 i? td 
Al y cr) sr? td ci 7 td 
AI Z ST? ST? td cr? td 




+A n 

+a 12 

+ Ai3 


Al x e x 

Xgd 


“ A J2 

+a 22 

+A 2 3 


Aly e y 



’ A 13 

- A 2 3 

+a 33 


_^Z e z _ 


where 


A a = +2(c 2 ij - s 2 r}) c 2 77 td = + c2t? (1 + c2r? td ) 


A12 = -2 ST? ST? td C7J td = - ST} s2T? td 
A13 = +2 CT} ST? td CT} td = + CT} s2T} td 

A 22 = +2(S 2 T} td - C 2 7} C 2 i? td ) = + (1 - c27? td ) - ~Y (1 + e2?}) (1 + c27?t d ) 
A 23 •= -2 ST} CT} c 2 t} td - - — s2t } (1 + c2T} td ) 

A 33 = -2(s 2 T? td - S 2 T} c 2 T} td ) = + I- ( 1 - c2t}) (1 + c2T} td ) - (1 - c2T} td ) 
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APPENDIX C 


COEFFICIENT EVALUATION 


The following definite integrals are involved in the development of the a^ 
coefficients of equation (7) or (8) 



s2i? td dr? td - 



s 2 7 td dlJ td = - -J O - c2 f>) 



(1 ± c27? td )dr? td 



(1 ± c27? td ) di? td - + -y- (2p ± s2p) 


With these integrals, we obtain ( 1 /2£2 has been absorbed in the Aj’s) 


&12 — 


0 0 
l J A 12 d7? td = -2sy J 


2 J A 12 = -2sr? J s 2 i 7 td d?? td - + srj (1 - c2p) 

-P “P 


a l3 “ 


1 


0 


0 


2 J A, 3 d? 7 td = +2cr ? J s 2 T 7 td di? td = - ct? (1 -c2p) 

-P -P 


a 23 




0 




0 


a 22 “ 


2j A 23 dr/ td = -s2 vj (1 + c2r? td ) dr? tc j = - -y s2t? (2p + s2p) 

-P -P 

0 0 0 

2 I A 22 dr? td = - (1 + c2t?)J (1 + c2r? td ) dr? td + 2 J (1 - c27? td ) dr? td 
-p -p -P 


= - -y ( 1 + c2i?) (2p + s2p) + (2 p - s2p) 
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0 

i 


di?td = + (1 + c2 n) / (1 + c2r) t( j) d7? td 


t 

i 


+ ~ (1 - c2r?) (2 p + s2p) - (2p - s2p) 


(1 - c2r? td ) dr? td 



APPENDIX D 


DESATURATION MANEUVER EFFICIENCIES 


The maneuver angles cannot be reached instantaneously.because of the limited 
vehicle angular velocity imposed by CMG momentum limitations. The angle profiles of 
Figure 3 are therefore used and the ratios between the desaturated momentum by the actual 
angle profile and the one desaturated by the ideal profile are expressed as efficiencies: 

E * = ( H dx) act / (»dx) id 
'V ■ ( H dy) ac t / ( H dy) id 
^ (Hdz) act / (H d z) id 

where the components are in CS Xp f . The actual desaturated momentum can be interpreted 

as the integration over infinitely small maneuvers with 100 percent efficiency, but with 
varying desaturation percentage. 

The ideal x desaturation momentum is 
( H dx) id - (1 - c2p) Ae' 


with 


Ae' = ST}Aey/(-A 2 ) - crjAe z /A 3 


[cf. equation (7)] . The actual x desaturation momentum is 



Ae' 

I [ 1 - c2p'] da 

* 0 
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where a is a dummy integration variable and 
p' = ap/(2Ae') 

Evaluation leads to 

Ae' 

(H dx ) acl = / o [-^(*-2f ? ) +B ^] d * 

= AsL f"s2p (l - + 

p [ l 2Ae / Ae'Jo 

I 

= ^(2sp-s2p) 

P 

and 

E = 2s P ~ s2 P (Dl) 

x p(a-c2p) 

The ideal y desaturation momentum is 

( H dy) id = 2s2 » € y 

with 

e y ' = C77e y /(-A 2 )-si7e z /A 3 

[see equation (9)] . The actual y desaturation momentum is 

Jy 

( H dy),c ' 2 s2( ' J '' , ' )d “ 
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with 


p = ap/( 2e y ') 
Evaluation leads to 



= (cp - c2p) 

P 


and 


F _ cp - c2p 

y ps2p 


(D2) 


Note that Ey is always larger than one (for the desaturation percentages considered) since 

the actual angle profile has less losses at the beginning and the end of the desaturation 
interval than the ideal profile (the torques are a function of c2r? td ). 

The ideal z desaturation momentum is 


( H dz) id = - (2p - s2p) e z ' 


with 


e z ' = si?e y /(-A 2 ) + C7?e z /A 3 
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[see equation (9)] . The actual z desaturation momentum is 


<H dz > 


act 


/ 1 2 (p - p ') - s2(p - p ')] da 
0 


where 


p' = ap/(2e z ') 


Evaluation yields 






(cp - c2p) 


1 e z' 


and 


E 


z 


|-p- — (cp - c2p) 

2 P 

2 p - s2p 


(D3) 
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APPENDIX E 


ORBITAL ELEVATION OF PRINCIPAL z p AXIS, NOMINAL ROTATION 
v zq ABOUT THE SUN LINE, AND GRAVITY 
GRADIENT TORQUE PHASE SHIFT n tm 

The sine and cosine functions of the orbital elevation j? of the principal z p axis 

are needed for equations ( 18 ), ( 19 ), and ( 20 ) and must be calculated from the known 
vehicle angles rj x and v z and the principal axes misalignment [K] . The transformation 

from orbital reference coordinates to principal reference coordinates is 


[c] = [K] [v z ] [tj x ] 


(El) 


which results in 


Cn * K n cv z - K 1;2 su z 

C,2 = (KuSV z + K 12 CV z )C 7 } x -K 13 S 7 ? x 

C13 = (K n Su z + K 12 co z )sr? x + Ki 3 Ci? x 
c 2 i = K 21 cu z -K 22 su z 

c 2 2 = (K 21 SU Z + K 2 2 Ci; z )CT? x - K 2 3Sr? x 
C 2 3 = (K 21 SU z + K22 cu z )s 1? x + K 23 C7 ? x 
C31 = K31 CV Z - K3 2 SU z 

C32 = (K31 SU Z + K 32 cu z )cr? x - K33 S1? X 
C33 = (K 31 SU Z + K 32 cu z )si? x + K33 cri x 

For the sine function of the elevation 7? of the principal z p axis we have (the attitude 
error is assumed zero) 
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SI? = -C32 = K33 S7? x - (K31 SU Z + K32 CU Z )ST? X 


(E2) 


The cosine function is 

ci? = V 1 - s 2 t?'SGN[K 31 su z + K 32 co z )si? x + K 33 ci? x ] (E3) 

There is no need to know the angle i? explicitly; the trigonometric functions are sufficient. 

The nominal i> z g rotation about the sun line puts the Xp axis into the orbital 
plane, or 

C12 = 0 

which yields 


K ji S» z g t CU z g k 1 3 ti?^ 


Kn K 13 ti? x - K12 VW + K12 2 ) - (K 13 ti? x ) 2 ' 
s «zg K u 2 + Kn 2 

_ K,2 K 13 tT? x + K„ /(K n 2 + K 12 2 ) - (K.3 ti? x ) 2 ' 
° Uzg K u 2 + K l2 2 


(E4) 


(E5) 


The phase shift 7 ? tm (midnight shift) of the cyclic gravity gradient torques (and 

consequently the shift of the resulting momenta) is the angle between the projection of the 
Xp axis into the orbital plane and the vehicle orbital x yo axis, which yields 


tr? tm = 


~Cl 3 

Cll 
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or 


’Jtm = 1 



-[(K,1 SU Z + K,2 cw z ).sjj x + K 13 C7J X ] 
K u cv z -K l2 sv 


} 


(E6) 
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APPENDIX F 


MANEUVER MOMENTUM PREDICTION AND ADAPTIVE 
ANGLE COMMAND LIMITING 


The maneuver commands of equations (21) and (22) are derived disregarding the 
fact that the CMG’s may not have sufficient momentum available to execute the maneuvers. 
A prediction of the available maneuver momentum must be made at the time the onboard 
digital computer calculates the maneuver commands, which must be scaled down, if neces- 
sary, to avoid the introduction of severe crosscoupling. 

The total CMG momentum H t to be predicted can be split into a part Hg that 

is unalterable (momentum caused by cyclic components, the average, and the ramp) and a 
part H man that can be changed (the momentum needed for the maneuvering itself and the 

momentum already desaturated by the maneuvers at the time in question): 


— t — g + man 


where X is a positive number (X = 1 indicates the unaltered case). The magnitude of this 
momentum is not allowed to exceed the saturation momentum H gat of the CMG’s: 

(Hg + XH man ) • (Hg + XH man ) - H sa £ 


or 


* = '(Hg ‘ —man) * V(H g * H man ) 2 + (H man . H man ) [H sat 2 - (H g . H g )]7 


(Hman * — man^ 


We can make the assumption that H sa £ > (Hg • Hg); i.e., without any maneuver the system 

is not saturated. Then only the plus sign in front of the square root results in a positive X 
(a negative value for X makes no sense; i.e., it indicates that a maneuver diametrically 
opposite to the one desired is necessary). When X > 1 , excess maneuver momentum is 
available and no scaling is necessary , whereas X < 1 requires scaling. 
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The number of time points where the maneuver momentum is checked must be 
minimized and the halfway points for each of the three maneuvers and the point at the 
very end of the maneuvers (end of the desaturation interval) have been selected as the check- 
points. This is illustrated in Figure F 1 which shows the orbital y and z components of 
the momenta, necessary for the maneuvering, added to the unalterable momentum, as a 
function of A different aspect of the same information is shown in Figure F2 where 

the z component is plotted as a function of the y component with 7 ?^ as a parameter. 

In this representation the saturation can be indicated as a circle about the origin in the 
y-z plane. The effect of an x momentum is also disregarded in Figures F 1 and F2 (for 
clarity), but the equations do not neglect the x component. In the following the X’s 
carry a subscript identifying the checkpoints with which they are associated. 

Saturation during the second maneuver is caused by a Ae only and can be detected 
at checkpoint 2. It requires a reduction of the Ac’s by X 2 if X 2 < 1. Checkpoint 2 is at 
the maximum momentum (for the second maneuver), if there is no orbital y momentum 
to be desaturated, and close to it for the case depicted in Figure FI (which has a large orbital 
y momentum). If saturation occurs at checkpoint 1 or 3, both the e’s and the Ac’s 
should be reduced proportionally by MIN(A 1 , X 3 ) such that saturation is just reached at 
the respective checkpoint. After the checkpoints, the desired momentum still can exceed 
the available momentum. Some momentum is always available right after the first maneuver 
and attaining the exact attitude is not critical; therefore, it makes no difference if the avail- 
able momentum is exceeded after its checkpoint. However, at the end of the third maneuver, 
the attitude error must be zero or experiment time will be lost. Figure F I shows that 
enough excess momentum is always available between the start of the third maneuver and 
checkpoint 3 to meet the time integral of momentum required for the desired attitude 
change. Increasing the commanded rate initially by a factor of p (which will be determined 
below) can be implemented by a fictitious reduction of the third maneuver interval by 
Ap. Figure F 1 shows that the problem can be linearized (which is slightly conservative) 
and the relationships are shown in Figure F3, where the available maneuver momentum is 
normalized with respect to the desired maneuver command. The cyclic momentum has been 
subtracted, resulting in a sloped saturation line, and the polarity has been reversed for con- 
venience. The momentum time integral requirement can be stated as (p is proportional to 
time) 

-j-M- y Pi (P-X 4 ) = -j- (FI) 

where it should be remembered that the ordinate is normalized by the desired maneuver 
momentum. Geometric relationships show that 

Pi = ~ (M-X 4 )/(X 3 -X 4 ) (F2) 


Elimination of p , from equation (F 1 ) yields 
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Figure F 1 . Checkpoints for maneuver momentum. 
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or 



Only the minus sign before the square root makes sense. To achieve a velocity increase by 
a factor of p by a shortening of the base by A p results in 

(f) f 


(F4) 

A plot of A p/p is shown in Figure F4. A conservative approximation that eliminates the 
need for a square root is 



~ 1 T 1 1 1 “ 1 1 1 I 

0 0.5 \ 1.0 

a 4 ► 

Figure F4. Exact Ap/(p/2) versus A 4 . 
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A plot of A p/p for the approximation is shown in Figure F5. All objectives of the sealing 
can be cast into three equations [see equations (18), (21), and (22)] 



(F6) 



MIN(1,\i,X 3 ) 



- MIN(1, A, , X 2 , A 3 ) 



(F7) 


Au z = -A 4 MIN( 1 , A 3 ) (sr)H dy + cr?H dz )/(27rcn x ) 


(F8) 
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For each checkpoint there is a need to identify what constitutes the nonscalable vector 
H n j and the scalable vector H se j. It is assumed that the components are in vehicle coordi- 
nates. In general [see equations (25), (G6), and (G7)] , 

— ni ~ — a + — b + ^nki— k + ^gyi— gy + ^gzi-gz + — dni (F9) 

where H^ni identifies the portion of the momentum already desaturated which is not 
affected by the particular Aj. Also 

Bsci = l'] li + Hda (F10) 

where ej signifies the change in angle during the maneuver in question and is the 

portion of the already desaturated momentum which is affected by the particular Aj. 

The following constants apply to the momentum ramp : 

K nki = 0.5 - 0.375 p/n 

K nk2 = 0-5 

K n k3 = 0.5 + 0.375 pin 
K n k4 = 0.5 + 0.5 p/w 

The following constants are associated with the cosine amplitude H g y (or H c ) 
and the sine amplitude H gz (or H s ) of the cyclic momenta [see equations (25), (G6), 
and (G7)] . Linearization was applied. 

Kgyi = - 1.5 p + 0.5 ir 

K gy2 = * 

K g y3 = - 1.5 p + 0.5 n 
K g y 4 — - 1 .8 p + 0.45 7r 

K g Z l = " 1 
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K gz2 = 0 


K gz3 = + 1 

Kg Z 4 — ■- 1 .7 p + 0.85 7T 

The desaturation momentum components must be in vehicle space, and the portion 
that will be desaturated by Av z after the desaturation interval (i.e., the orbital z momen- 
tum) must be subtracted: 


(F 11) 


The portion of the momentum that must be considered unscalable is 


-dvx 


^dx 


r “i 

0 

^dvy 

= IK] T / 

Hdy 

- (srjH dy + cr?H dz ) 

ST? 

^dvz 

( 

H dz 


CT? 


H dn i = [0,0, 0] T 


H dn2 = 0.5MIN(l,X 1 ,X 3 )[0,H dvy ,H dvz ] T 


H d n3 


— dn4 


[0, 0, 0] T 

MIN( 1 , X, , X 2 , X 3 ) H dvx 

MINCUX^Xa)^^ 

MIN(l,X 1 ,X 3 )H dvz 


The portion of the desaturation momentum that can be scaled is 


H dsl = [0, 0, 0] T 
Hds2 = l0.5H dvx ,0,0] T 
— ds3 _ ^dvx> ^dvy> ^ dvz ]^ 

H d8 4 = 10,0,0] T 
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The attitude changes during the maneuvers are 


—1 ^ e xl’ e yl’ e zO 
e_2 ~ [-Ae x , -Ae y , -Ae z ] T 

e 3 ~ f' e x2’ ~ e y2’ ’ ( e z2 " Au z^ ^ 

where equations (18), (19), and (20) supply the e’s and the Ac’s (i.e., the angles before 
scaling) but 

—4 = [- e x2’- e y2’^ e z2- Au z)] T 

where equations (F7) and (F8) supply the e’s and Av z (i.e., the angles after scaling). 
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APPENDIX G 


CYCLIC MOMENTUM COMPONENTS 


The cyclic momentum components can be extracted from the samples [equation 
(23)] . However, in some cases one or more samples are invalid or missing, and an alternate 
way of determining the cyclic momentum components must be available. 

The intended use is not critical, and it can be assumed that the large moments of 
inertia (about the yp and Zp axes) are identical. This assumption allows a rotation about 

the Xp axis without a change in the gravity gradient torque vector. It is then convenient 
to develop the torque with the Xp axis out of the orbital plane by an angle u ze and with 
the Zp axis in the orbital plane: 


T g = 3 £2 2 AIr x 


with 



cu ze 

su ze 

0 

su ze 

CtJ ze 

0 

0 

0 

1 


CT7 t d 0 S7J td 


" 0 " 

i 

0 1 0 


0 

-S77 td 0 cij td 


-1 


-cv 


ze 


su 


ze 


S7 ?td 

S7 ?td 


' c,J td 


A I ‘ Ox ‘ I z > iy " ^x 
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or 

r o 


Xg ■+ ~Y O 2 AIcv zq 


s2 7Jtd 
[-su ze (l - c 2 i 7 td ) 


(Gl) 


Since we are only interested in the cyclic components, the bias in the z component can be 
ignored: 


0 




£2 2 Alcu ze 


-s2r ?td 


(G2) 


su ze c2r »td 


Integration yields 



-f- “ Alcu ze 


0 ' 
c2 ^td 


su ze s2? ?td 


(G3) 


where the average over one orbit has been set to zero. To be useful the cyclic momentum 
components should be in vehicle components: 


— gcv 


- £2 Alcu ze [u 7 ] [?? Y ] 


tm J 


"ze J 


0 

c2r ?td 
su ze s2p td 


(G4) 


where r? tm is the angle about orbital y between the projection of the Xp axis into the 
orbital plane and the z y axis into the orbital plane. Evaluation of equation (G4) yields 
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— gvc — gy C ^ TJ td + — gz s ^td 


(G5) 


with 



3 

4 


£2 Alc» ze 


-CV Z Cl7^ m Sl» ze + SU Z (cr? x CU ze + ST? X STJ j-jjj SJ> ze ) 

su z C7J tm su ze + cu z ( CT »x cw ze + S7 ?x S7 »tm sw ze> 
-S7j x cu ze + ct? x si? tm su ze 


(G6) 


Hgz = Y flAls2u ze 


cu z ST »tm + su z ST tx CT »tm 
-s« z srj tm + cv z s Vx cv tm 

CT ?x C7? tm 


(G7) 


Only the terms s2i7 td and c2ij t( j must be updated frequently ; the rest can be updated once 

per orbit. If i> ze and r?j m can be considered small, equations (G6) and (G7) can be sim- 
plified to 


Hgy = 4 nal 


-° v z w ze + sv z CJ ?x 
so z u ze + cv z CT »x 

-sjj x 


(G8) 


Hg Z = T 


SU Z STJ X 
co z stj x 

CTJ X 


'ze 


An estimate must be made on u, 


ze - 


u ze ( u z ' u zg ) CT? x 


(G9) 
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where u zg is the z rotation needed to put the x p axis into the orbital plane (see Appendix 
E for its derivation). Further simplification is possible if o ze can be considered negligible: 




” sv z CT ?x ~ 

CU Z C7? x 



(G10) 


0 



0 

0 


(G 11) 


Hgy of equation (GIG) is the resolution of a momentum along the orbital y axis into vehicle 
components. 
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APPENDIX H 


SAMPLED DATA SYSTEM 


The following equations are needed to treat the desaturation method as a sampled 
data system: 


H n = + 


(HI) 


EH n 2H n _j + K n H n + K n _j H n _j 


(H2) 


^dn ^ ^n ' 


(H3) 


[see equations (28) and (29)] where H is the total accumulated momentum and H a = H. 
% is the momentum accumulation per orbit. 

Application of the z transform results in 


h = z 1 h + h k + z 1 hj 


2h = z _1 2h + K n h + z' 1 K n _jh 


h d = 2h-h k 


(H4) 

(H5) 

(H6) 


or 


h = (z- l) 2 

h k z 2 -(2 + K n )z + (l-K n .j) 

The characteristic equation 


(H7) 


z 2 -(2 + K n )z + (l-K n . 1 ) 
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results in 


z h2 = 0.5 [(2 + K n ) ± V(2 + K n ) 2 -4(1 -K n .,)] 


(H8) 


K n is selected to be -1, which yields 


z j 2 = 0.5 ± VK n .! - 0.75 


This results in an oscillatory system for 


K n .j < 0.75 


and the equivalent system 1 to equations (27) for 


K 


n-1 “ 


1 


since the z = 1 is cancelled by the numerator in H7. 


1. Only true as long as no nonlinearities are encountered. 
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